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THE  DESIGN  AND  ANALYSIS  OF  A  UNIQUE 
BROADBAND  UNDERWATER  ACOUSTIC  SOURCE 


INTRODUCTION 

Underwater  electroacoustic  calibration  measurements  taken  in  open  bodies  of  water 
are  typically  made  at  relatively  shallow  depths.  To  minimize  the  multipath  effects  of  reflec¬ 
tions  from  the  surface  and  bottom,  it  is  often  desirable  to  use  a  calibration  source  that  dis¬ 
criminates  against  those  directions  over  the  frequency  range  of  interest;  that  is,  a  source  that 
essentially  confines  the  radiation  within  a  relatively  narrow  '  am  in  the  vertical  plane.  Cali¬ 
bration  measurements  can  also  be  simplified  if  the  output  of  the  source  is  constant  over  the 
frequency  range  in  which  it  will  be  used.  However,  this  report  shows  that  these  two  require¬ 
ments  are,  in  general,  contradictory. 


REQUIREMENTS 

Directivity  is  typically  the  most  difficult  parameter  to  control  over  a  wide  frequency 
range.  The  intended  application  for  this  transducer  requires  a  toroidal  beam  pattern;  that  is, 
omnidirectional  in  the  horizontal  plane  with  a  main  beamwidth  between  30  to  60°  (at  the 
3-dB  down  points)  in  the  vertical  plane.  The  sidelobes  are  required  to  be  a  minimum  of 
15  dB  below  the  maximum  level  of  the  main  beam.  The  directivity  requirements  are  sum¬ 
marized  in  Fig.  1. 

The  directivity  requirements  of  the  subject  transducer  determine  its  basic  geomet¬ 
rical  configuration.  The  requirements  for  omnidirectionality  in  one  plane  and  a  specified 
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beamwidth  in  a  perpendicular  plane  imply  the  geometry  of  a  linear  array  with  the  longi¬ 
tudinal  array  axis  perpendicular  to  the  omnidirectional  plane.  Also  implied  by  the  omni¬ 
directionality  in  the  horizontal  plane  is  the  cylindrical  symmetry  of  the  individual  radiating 
elements.  Therefore,  based  upon  the  directivity  requirements,  it  can  be  assumed  that  the 
subject  transducer  will  be  in  the  basic  form  of  a  linear  array  of  radially  poled  cylindrical 
piezoelectric  ceramic  shells. 

The  system  incorporating  the  subject  transducer  requires  an  output  sound  pressure 
level  (SPL)  of  190  dB  (referenced  to  1  pPa  measured  at  1  m)  over  the  frequency  range 
from  10  to  40  kHz.  It  is  required  that  the  transmitting  voltage  response  (TVR)  be  constant 
(±1.5  dB)  over  the  same  frequency  range.  The  typical  TVR  of  a  piezoelectric  ceramic  cylin¬ 
drical  shell  with  the  ends  and  inner  surface  acoustically  shielded  (or  a  linear  array  of  such 
elements)  is  not  flat,  but  instead  has  a  +12  dB  per  octave  slope  at  frequencies  well  below  the 
frequency  of  the  first  radial  resonance.  Obviously  then,  a  simple  linear  array  of  cylindrical 
elements  will  not  meet  the  TVR  requirements. 

The  requirement  for  the  flat  response  demands  that  the  distribution  of  the  volume 
velocity  generated  by  the  array  be  such  that  the  sum  of  the  pressures  in  the  far-field  is  a 
constant.  This,  coupled  with  the  directivity  requirements  that  the  array  length  be  a  function 
of  frequency,  implies  the  biconical  configuration  shown  in  Fig.  2.  The  elements  are  arranged 
symmetrically  with  respect  to  the  array  center  according  to  the  frequency  of  their  radial 
resonance;  that  is,  the  largest,  lowest-frequency  elements  are  at  the  ends  of  the  array  with 
the  element  size  decreasing  and  the  resonance  frequency  increasing  symmetrically  toward 
the  array  center.  When  the  elements  are  driven  electrically  in  parallel,  the  effective  array 
length  is  a  function  of  frequency  with  the  array,  in  genera),  becoming  shorter  with  increas¬ 
ing  frequency. 

In  the  succeeding  sections,  this  report  documents  the  development  of  a  model  of  a 
generalized  linear  array  of  cylindrical  acoustic  radiators;  the  design,  fabrication,  and  evalua¬ 
tion  of  a  prototype  of  the  subject  transducer;  and  the  comparison  of  the  measured  results 
with  those  predicted  by  the  model. 
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Fig.  2  —  The  biconical  array  configuration  implied  by 
the  TVR  and  directivity  requirements 
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ANALYSIS  OF  THE  SUBJECT  TRANSDUCER 

In  terms  of  what  is  already  known  about  linear  arrays,  the  desired  transducer  will  be  an 
unequally  spaced,  amplitude  and  phase-shaded  linear  array  [1-6] .  The  elements  will  be  un¬ 
equally  spaced  primarily  because  of  the  requirement  to  tailor  the  volume  velocity  distribu¬ 
tion  to  produce  the  flat  TVR;  that  is,  different  element  lengths  are  required  for  each  ele¬ 
ment  pair  based  upon  their  radial  resonance  frequency.  Since  the  elements  are  connected 
electrically  in  parallel,  the  same  voltage  amplitude  is  present  across  each  of  the  elements  in¬ 
dependent  of  the  frequency.  The  array  will,  however,  be  amplitude  shaded  for  two  reasons. 
First,  at  any  one  frequency  there  is  a  relative  difference  in  element  amplitude  due  to  their 
different  frequencies  of  radial  resonance.  Second,  the  “adjusted”  volume  velocity  distribu¬ 
tion  necessary  for  the  required  TVR  is,  of  course,  a  form  of  amplitude  shading.  It  should 
be  obvious  that  the  amplitude  shading  required  by  the  transmitting  response  and  the  direc¬ 
tivity  is  not  necessarily  the  same.  In  fact,  the  two  requirements  are  for  the  most  part  contra¬ 
dictory.  The  array  will  be  phase  shaded  because,  for  a  given  frequency,  the  response  curves 
of  the  individual  element  pairs  are  at  different  points  relative  to  their  resonance  frequencies; 
in  other  words,  there  is  a  phase  difference  between  the  surface  velocities  of  the  element 
pairs. 

The  biconical  configuration  is  essential  to  the  directivity  requirement  in  that  the  array 
will  appear  to  effectively  shorten  as  a  function  of  increasing  frequency.  The  key,  of  course, 
is  controlling  the  response  of  each  element  pair  at  frequencies  above  resonance.  One  simpli¬ 
fying  alternative  would  be  to  electrically  series  tune  each  element  pair  to  its  frequency  of 
radial  resonance.  However,  to  do  so  would  severely  limit  the  bandwidth  obtained  from  each 
pair  and  would  result  in  a  requirement  for  many  more  closely  spaced  elements,  both  physi¬ 
cally  and  in  frequency.  In  such  a  configuration  the  very  small-diameter,  thin-walled  cylin¬ 
ders  required  for  the  highest  frequencies  could  not  produce  the  required  SPL. 

The  biconical  configuration  also  results  in  a  reverse  shading  or  negative  tapering  of  the 
array.  For  example,  at  the  resonance  frequency  of  the  largest,  lowest-frequency  elements 
the  amplitude  is  highest  at  the  ends  of  the  array  and  decreases  toward  the  center.  As  might 
be  expected,  this  condition  leads  to  a  narrower  main  beam  and  higher  sidelobe  levels  1 7  ] . 
This  may  be  partially  overcome  by  judiciously  choosing  the  eenter-to-center  spacing  of  the 
elements. 

In  summary,  the  primary  design  parameters  (all  of  which  are  interrelated)  are  the  effec¬ 
tive  array  length  (as  a  function  of  frequency);  the  eenter-to-center  spacing  of  the  element 
pairs;  the  radiating  area  of  the  individual  element  pairs;  and,  of  course,  the  relative  positions 
(in  frequency)  of  the  radial  resonance  frequencies  of  the  elements. 


ANALYSIS 

For  the  sake  of  simplicity,  consider  first  the  case  of  the  far-field  pressure  generated 
by  an  array  of  two  point  sources  radiating  in  phase.  The  generated  pressure  will  be  of  the 
form 
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iO  if) 

l+v 


(1) 


3 


YOUNG 


where  A 1  and  A2  are  the  respective  pressure  amplitudes  and  0  j  and  02  represent  the  phase 
of  the  two  points  with  respect  to  the  geometrical  center  of  the  array.  Looking  only  at  the 
real  pressure  amplitude, 

Pj  2(0)  =  £(-Ai  cos  81  +  A2  cos  d2)z  +  {A1  sin  61  +  A2  sin  02)2J  ^  >  (2) 

or 

PlM0)=Al  +A2  +2A1A2  COS(01  ~  °2>’  (3> 

where  the  argument  of  the  cosine  term  01-  02  represents  the  total  phase  difference  between 
the  two  points.  From  the  geometry  of  Fig.  3,  it  is  obvious  that  Q1  =  kx  and  02  =  ~kx,  where 
k  =  2rr/X  and  X  is  the  acoustic  wavelength  in  water.  Thus  it  is  a  simple  matter  to  solve  for  x 
to  yield  =  fed/2  and  02  =  (-fed/2)  sin  0. 

If  we  now  consider  the  two  sources  to  be  cylindrical  elements  of  finite  but  unequal 
radii,  the  phase  relationship  between  the  two  is  shown  in  Fig.  4.  Now  the  total  phase  differ¬ 
ence  between  the  sources  due  to  their  separation  in  space  is  k(x  -  r),  where  x  and  r  may  be 
found  from  Fig.  4  to  be  d  sin  0  and  ( a2  -  )  cos  6,  respectively.  If  the  driving  function  (the 

electrical  signal)  is  applied  to  the  two  elements  in  parallel,  another  phase  difference  becomes 
apparent. 

Since  the  elements  do  not  have  the  same  mean  radii,  their  radial  resonances  will  occur 
at  different  frequencies  and  their  surface  velocities  will  therefore  differ  in  phase.  If  the 
phase  of  the  surface  velocity  with  respect  to  the  driving  function  is  <t>1  and  <p2,  respectively, 
Eq.  (3)  for  the  square  of  the  generated  far-field  pressure  becomes 

Pf  2(0)  =  A  j  +  A\  +  2.A  ji42  cos  |  8j  -  02  -  fe  j^d  sin  0  -  (a2  "  ai )  cos  ^  \  ■  (4) 


Fig.  3  —  An  array  of  two  point  sources  separated  by  a  distance  d 
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Fig.  4  --  An  array  of  two  cylindrical  sources  of  unequal  radii 
separated  by  a  distance  d 


The  square  of  the  on-axis  (6  =  0)  pressure  reduces  to 

Pf  2(0)  =  Af  +  A%  +  2A1A2  cos  -  <p2  +  k(a2  -  Oj )  J  ,  (5) 

and  the  normalized  directivity  pattern  becomes 


•*\Z (0) 


P1  2(0)  (  A\  +  a\  +  2A}A2  cos  (0j  -  <t>2  -  (Aj  2)  yu 

^1,2^  j  +  a\  +  2^1, i42  cos  £</>,  -  <t>2  +  fe(a2  -  dj  )J  \ 


(6) 


where  t//  j  2  =  fe  sin  0  -  (a2  "  ai )  cos  • 

The  elements  are  now  of  finite  dimensions  and  can  no  longer  be  considered  as  point 
sources;  that  is,  there  is  a  nonunity  directivity  factor  associated  with  each  element.  The 
pressure  amplitudes  of  the  individual  elements  as  a  function  of  0  now  can  be  represented  by 
flj  [sin  (fetj/2'sin  f?)/(/eC1  /2-sin  0)]  andB2=^2  [sin  (fet2/2*sin  0)/(fel!2/2*sin  (1)] 
where  the  right-hand  term  in  each  expression  is  the  normalized  directivity  pattern  for  a  con¬ 
tinuous  line  source  of  length  Cj  and  C2>  respectively.  This  can  be  shown  to  be  simply  an 
application  of  the  product  theorem  [8] ,  which  (as  applied  to  this  case)  states  that  if  the 
point  sources  in  the  two-element  array  are  replaced  by  finite  sources,  then  the  resulting 
directivity  pattern  for  the  array  will  be  the  product  of  the  directivity  patterns  of  the  finite 
elements  and  the  two-element  point  source  array.  To  illustrate,  Eq.  (6)  may  be  rewritten  as 


•'l.2<0> 


\  B\  +  fl|  +  2BlB2cos(<t>^  -02-  *12)  |  1/2 

j  A  J+  A I  +  2A  j/t2  cos  jVj  -  <t>2  +  al  \ 


(7) 
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and  if  the  elements  are  of  equal  size,  Eq.  (7)  then  reduces  to 


*^i  2 


■«  i  r 

sin  —  sin  0  sin  kd  sin  0 

L  2  J  L  J 


■fee 

'fed 

—  sin  0 

L 2 

2  sin 

• —  sin  d 
.  2 

(8) 


which  is  the  product  of  the  directivity  patterns  of  a  continuous  line  element  of  length 
with  that  of  an  array  of  two  point  sources  separated  by  a  distance  d. 

This  overall  approach  may  be  generalized  to  include  N  arbitrarily  dimensioned  and 
spaced  finite  elements,  and  Eq.  (4)  may  be  written  as 


C°S  (0„  -  <P„,  ~ 


li/  ) 

Y  n.m  '  ’ 


(9) 


where  Bn  =  An  sin  (feC„/2  •  sin  d)/(kin/2  •  sin  0), 

Bm  =  Am  sin  /2  '  sin  d)Kk^m  /2  sin  #). 

'/'n.m  =  k  ldn,m  sin  0  ~  (am  an  )  cos  0]  .  and 

dn  m  (or  dmj7 )  is  the  spacing  between  the  nth  and  wth  elements. 

In  generalized  form,  the  square  of  the  on-axis  pressure  is  written  as 

Pn(°)  =  ^  J2  AnAm  COS  fy  Ki  +  -  an>]  - 

n«l  m  = 1  L  J 

and  the  directivity  pattern  becomes 


cos  (<b  -  <j>  -  \b  ) 

n  'Yn  Ym  Y n.m  ’ 


AnAm  COS[^  ~  +  ~  °n  > J 


(10) 


(ID 


Before  these  expressions  can  be  used  to  compute  the  TVR  or  the  directivity  patterns  in 
the  vertical  plane,  the  pressure  amplitudes  (the  A  coefficients)  for  the  individual  cylindrical 
elements  must  be  determined.  In  general,  the  far-field  pressure  from  each  element  can  be 
represented  by  [9] 


P  = 


\  u2rarqPqc  { 
)  4nr2  ^ 


1/2 


(12) 
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where  V  is  the  root-mean-square  (rms)  volume  velocity  of  the  radiator,  R A  is  the  acoustic 
resistance  acting  on  the  radiator,  Rq  is  the  directivity  factor,  p0c  is  the  characteristic  im¬ 
pedance  of  the  medium,  and  r  is  the  distance  to  the  measurement  point.  The  diffraction 
constant  for  such  a  source  may  be  written  as  [9] 


J  RAR6*v  11/2 

I  k2p0c  i 


(13) 


which  may  be  solved  for  Rq  and  substituted  into  Eq.  (12).  With  this  substitution,  the 
expression  for  the  pressure  becomes 


UDkpQc 

P=— -  .  (14) 

4  irr 

Using  the  definition  of  the  volume  velocity  (U  *  u  •  S)  and  the  surface  area  of  the 
cylinder  (S  =  27T aft),  the  pressure  may  be  expressed  as 

upQcDkaQ 

F - iT-1  ,16) 

where  u  is  the  rms  velocity  of  the  radiating  surface,  a  is  the  radius  of  the  cylinder,  and  C  is 
the  cylinder’s  length.  Now  it  is  only  necessary  to  find  an  expression  for  the  surface  velocity 
in  order  to  compute  the  far-field  pressure  produced  by  each  element. 

The  radial  velocity  of  the  cylindrical  elements  will  be  of  the  form 

V  =  u  cos  (cot  -  0)  ,  (16) 

where  the  velocity  amplitude  u  is  the  ratio  of  the  rms  force  to  the  mechanical  impedance 
of  the  ceramic  cylinder  and  <p  is  the  phase  difference  between  the  driving  function  and 
velocity.  The  mechanical  impedance  and  phase  angle  are  defined  conventionally  as 
Z,n  =  Rm  +  lRm  80111  ^  =  tan  1  lRm  )  where  Rm  and  Xm  are  the  mechanical  resis¬ 
tance  and  reactance,  respectively.  By  using  the  definition  of  the  mechanical  quality 
factor  Q  and  the  fact  that 


X  =  coM  -  — - 
m  wC. 


the  magnitude  of  the  impedance  is  found  to  be 


(17) 
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where  M  and  Cm  are  the  mass  and  mechanical  compliance,  gj  is  the  frequency  in  radians  per 
second,  toQ  is  the  radial  resonance  frequency  in  radians  per  second,  and  Q  =  1/c oQRm  Cm  . 
Similarly,  the  phase  angle  can  be  found  to  be 


The  rms  force  produced  by  the  piezoelectric  ceramic  is  given  by 

27rCd31  V 


(18) 


(19) 


d31  is  the  piezoelectric  strain  constant,  where  the  subscripts  denote  the  axes  of  the  applied 
electric  field  and  induced  strain  respectively,  S^j  is  the  reciprocal  elastic  modulus  at  constant 
electric  field  (Sfj  is  the  reciprocal  of  Young’s  modulus),  and  V  is  the  rms  value  of  the  ap¬ 
plied  voltage.  If  we  now  take  the  ratio  of  the  magnitude  of  the  force  to  the  magnitude  of  the 
mechanical  impedance,  we  get 


2JrMslV«Cm<? 


u  = 


oE 

*11  | 


<  "2  .  .2 


1/2 


(20) 


for  the  magnitude  of  the  velocity.  Substituting  Cm  =  aSf,  2n\.t  1 10)  end  H m  -  p0O27raV ) 
along  with  the  above  expression  and  the  expression  for  tan  0  into  Kq  1 1 5 1  and  performing 
the  algebraic  manipulation,  we  obtain 


d31  V(p/Sfj  )1/2IM'uho0  cos  0 


for  the  pressure  generated  by  the  individual  elements  (where  o  is  the  density  of  the  piezo¬ 
electric  ceramic).  Equation  (21)  may  be  substituted  into  Eq.  (9)  for  the  .4  coefficients 
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and  both  sides  of  the  result  divided  by  V2  to  yield  the  far-field  pressure  per  volt  from  the 
generalized  array.  The  result  is 


'N(0)|  COdg^p/Sfj) 


E  vl/2 


a;  *  Kn  sin  un  Km  sin  um  | 

Z  Z  T~  — -  COS  "  K  -  *«.*  )  .  (22) 

n  =  1  m  =  1  7 


where  Kn  =  DnanVn  cos  <fin,Km  =  Dmam£mcjm  cos<j>m,un  =  (kUn  /2)  sin  G,  and 
um  =  (feCm  12)  sin  0.  The  on-axis  sum  of  the  pressures  per  volt  (the  TVR)  reduces  to 


PN( 0)  udn(p!S^)112  ( 


^  KnKm  cos  [K  "  K  +  fe(“m  -  aj] 


and  the  directivity  pattern  of  the  N -element  line  array  becomes 


N  K„  sin  u„  K„  sin  u„ 
ir^  n  n  m  m 


•*V0) 


n  n  m  m  J 

Z  Zt~~  ;  coswn~0m 

n  =  1  m  =1  n  m  \ 

|  N  N  ( 

Z  Z  COS  +  *<am  "  V  |  ) 

n=  1  m  =  1  >-  ' 


The  cylindrical  elements  have  been  assumed  to  vibrate  only  in  the  radial  mode,  and  the 
mutual  effects  between  radiating  elements  have  been  neglected,  or  at  least  to  this  point 
lumped  into  one  unknown  parameter— the  diffraction  constant. 

If  the  mutual  coupling  effects  between  elements  are  ignored,  the  precise  diffraction 
constant  for  the  biconical  configuration  is  still  unknown  and  would  be  extremely  difficult 
to  determine.  However,  expressions  do  exist  for  the  diffraction  constants  of  a  thin  cylin¬ 
drical  ring  and  a  long  cylinder.  By  direct  integration,  Henriquez  [11]  found  the  diffraction 
constants  for  these  two  configurations  to  be 

Dring='/0<*a>  (25) 


^cylinder  nka 


Jj2(fea)  +  N2(fea)J 
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where  J0  and  Jj  are  the  zeroth  and  first  order  Bessel  functions  and  jVj  is  the  first  order 
Neumann  function.  Considered  individually,  the  array  elements  are  probably  best  described 
as  thin  rings.  In  the  array,  however,  adjacent  elements  act  as  a  tapered  cylindrical  baffle  (at 
least  geometrically)  and  the  long  cylinder  configuration  is  perhaps  more  descriptive.  Equa¬ 
tion  (26)  will  be  used  as  an  approximation  to  the  diffraction  constants  of  the  array  elements 
in  Eqs.  (23)  and  (24). 

A  program  was  written  to  carry  out  the  calculations  indicated  in  Eqs.  (23)  and  (24)  on 
a  PDP-11/45  digital  computer.  Data  input  is  required  to  specify  the  piezoelectric  material 
type  (the  piezoelectric  constant,  the  reciprocal  elastic  modulus,  and  the  material  density), 
the  element  dimensions  (outside  diameter,  wall  thickness,  and  length),  and  the  center-to- 
center  spacing  of  the  elements  in  the  array.  There  is  the  option  of  selecting  the  program  out¬ 
put  as  either  the  TVR  over  a  specified  frequency  range  or  the  directivity  patterns  at  speci¬ 
fied  discrete  frequencies.  The  data  may  be  displayed  in  tabulated  form  on  a  terminal  or 
plotted  on  an  X-Y  plotter. 


THE  PROTOTYPE  ARRAY 

To  provide  a  relatively  simple  first  test  of  the  model,  a  prototype  array  was  designed 
to  operate  over  approximately  half  of  the  total  required  bandwidth.  The  array  contains 
four  elements  (two  element  pairs)  and  operates  from  15  to  30  kHz.  An  iterative  procedure 
using  the  computer  program  to  compute  the  TVR  was  utilized  to  determine  the  element 
sizes.  The  element  diameters  and  wall  thicknesses  were  varied  within  the  constraints  of  avail 
able  ceramic  sizes  and  the  lengths  varied  to  obtain  the  desired  flatness  of  response.  While 
maintaining  the  same  ratio  of  element  lengths  and  remaining  within  the  constraints  of  the 
required  SPL  and  electroacoustic  efficiency,  the  elements  were  then  shortened  in  order  to 
minimize  the  element  spacing  in  the  array.  Using  this  procedure,  the  array  shown  in  Fig.  5 
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was  derived.  The  center-to-center  spacing  of  the  elements  is  not  simply  one  half  of  the  sum 
of  the  lengths  of  adjacent  elements,  because  of  the  space  occupied  by  mounting  hardware. 

The  generalized  model  describes  an  array  where  the  elements  are  radially  vibrating  cy¬ 
lindrical  radiators  free  from  any  mechanical  or  acoustical  coupling  effects.  We  have  already 
chosen  to  ignore  any  mutual  acoustical  effects  to  simplify  the  analysis,  but  care  must  be 
taken  in  the  design  to  assure  that  the  elements  are  not  mechanically  coupled  through  the 
mounting  structure.  To  mechanically  decouple  each  element  from  its  mounting,  glass-loaded 
polycarbonate  end  rings  are  fastened  to  each  end  of  the  ceramic  cylinders  with  a  thin  film 
of  relatively  compliant  potting  compound.  Since  the  elements  are  to  be  air  backed,  the  end 
rings  seal  against  a  center  mounting  spindle  by  using  elastomer  O-ring  seals.  The  spindles  are 
in  turn  mounted  to  a  single  central  shaft  in  the  transducer  by  a  cast  ring  of  compliant  pot¬ 
ting  compound.  A  single  element  and  its  mounting  configuration  are  shown  in  Fig.  6. 


-  RETAINING 
PLATE 

"0"  RING 


PIEZOELECTRIC 
-  CERAMIC 
CYLINDER 


POTTING 

COMPOUND 


Fig.  6  -  An  exploded  view  of  a  single  prototype  element  and 
its  mounting  structure 
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Cylindrical  stainless  steel  end  plates  are  fastened  to  each  end  of  the  central  mounting 
shaft,  and  the  transducer  is  sealed  in  an  elastomer  cylindrical  boot  and  filled  with  castor  oil. 
The  complete  transducer  is  shown  in  Fig.  7. 


COMPARISON  OF  RESULTS 

The  acoustical  characteristics  of  the  prototype  array  were  measured  in  the  frequency 
range  from  15  to  30  kHz  and  to  hydrostatic  pressures  of  4  MPa.  Measurements  made  in¬ 
cluded  TVR,  directivity  patterns  in  the  horizontal  (XY)  and  vertical  (XZ)  planes,  and  the 
linearity  of  output  SPL  as  a  function  of  driving  voltage.  All  of  the  measurements  show  very 
little  change  (less  than  0.5  dB)  as  a  function  of  hydrostatic  pressure,  and  the  horizontal 
directivity  patterns  are  omnidirectional  (±1.0  dB).  The  transducer  is  linear  as  a  function  of 
driving  voltage  within  ±0.5  dB. 

The  array  easily  produces  the  required  on-axis  SPL.  At  an  output  SPL  of  190  dB 
(reference  to  1  //Pa  and  measured  at  a  distance  of  1  m)  the  electric  field  on  the  smallest, 
thinnest-wall  ceramic  is  1000  V/cm,  or  approximately  half  of  the  voltage  it  can  safely 
withstand. 

The  TVR  as  predicted  by  Eq.  (23)  is  compared  to  the  measured  response  in  Fig.  8. 

The  two  response  curves  agree  to  within  1.5  dB  or  less  and  are  flat  (±1.5  dB)  over  the  design 
frequency  range  of  15  to  30  kHz.  The  close  agreement  between  the  curves  indicates  that  the 


Fig.  7  —  The  assembled  prototype  transducer 
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Fig.  8  —  Comparison  of  the  measured  TVR  of  the  prototype  array  with 
that  predicted  by  the  model 


error  introduced  by  neglecting  mutual  radiation  effects  between  the  elements  is  small,  at 
least  for  the  case  of  the  prototype  transducer. 

The  predicted  (Eq.  24)  vertical  directivity  patterns  are  compared  to  the  measured 
patterns  at  15,  20,  and  30  kHz  in  Figs.  9  (a),  (b),  and  (c)  respectively.  Only  half  of  the  full 
vertical  patterns  is  shown  in  the  figures  because  the  other  half  is  simply  its  mirror  image. 

Some  general  conclusions  can  be  reached  about  the  measured  patterns  in  terms  of  the  speci¬ 
fications;  the  width  of  the  main  beam  at  the  3-dB  down  points  is  marginally  acceptable  but 
the  sidelobe  levels  are  too  high,  particularly  at  the  lower  frequencies.  Some  general  observa¬ 
tions  may  also  be  made  about  the  accuracy  of  the  model:  while  the  width  of  the  main  beam  is 
reasonably  well  predicted,  the  fine  structure  of  the  pattern  (the  nulls  and  sidelobe  levels)  is 
not.  The  prediction  of  the  width  of  the  main  beam  is  reasonable,  probably  because  the 
diffraction  constant  used  (Eq.  26)  and  the  normalized  directivity  pattern  assumed  for  the 
individual  elements  in  the  model  are  accurate  for  the  narrow  angular  confines  of  the  main 
beam.  As  the  observation  point  becomes  further  off  the  axis  of  the  main  beam,  the  unac¬ 
counted  for  effects  form  the  adjacent  elements  and  the  housing  end  plates  alter  the  effec¬ 
tive  diffraction  constants  and  the  normalized  directivity  patterns  of  the  individual  elements. 
This  is  probably  best  illustrated  by  examining  the  measured  directivity  pattern  at  15  kHz. 

The  sidelobes  of  the  pattern  look  very  similar  to  what  one  would  expect  from  the  larger 
two  end  elements  spaced  0.8  X  apart;  that  is,  the  effect  of  the  center  elements  is  less  than 
predicted  by  the  model. 


CONCLUSIONS 

The  directivity  characteristics  of  the  prototype  array  can  be  brought  to  within  the 
specifications  over  the  15-  to  30-kHz  frequency  range  very  easily;  that  is,  the  array  can 
simply  be  shortened  to  an  effective  length  of  approximately  X/2.  This  is  a  feasible  approach 
in  the  case  of  the  prototype  because  we  can  reduce  the  ceramic  volume  by  almost  half  and 
still  maintain  a  safe  driving  voltage  level.  If  the  lengths  of  the  element  pairs  are  scaled  from 
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330°  6°  30° 


Fig.  9(a)  —  Comparison  of 
the  measured  vertical  (XZ) 
directivity  pattern  at  15 
kHz  with  that  predicted 
by  the  model 


Fig.  9(b)  —  Comparison  of 
the  measured  vertical  (XZ) 
directivity  pattern  at  20 
kHz  with  that  predicted 
by  the  model 
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Fig.  9(c)  —  Comparison  of 
the  measured  vertical  (XZ) 
directivity  pattern  at  30 
kHz  with  that  predicted 
by  the  model 
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the  present  configuration,  the  flatness  of  the  TVR  will  be  maintained  but  the  output  level 
per  volt  will,  of  course,  be  lower.  Reducing  the  volume  of  piezoelectric  ceramic  will  also 
lower  the  electroacoustic  efficiency  although  some  of  the  loss  may  be  regained  by  parallel 
tuning  the  array  with  the  appropriate  inductor. 

Model  predictions  indicate  that  an  array  designed  for  the  full  required  frequency  range 
of  10  to  40  kHz  will  require  eight  elements  (four  element  pairs).  The  safe  driving  voltage 
limit  for  such  an  array  is  determined  by  the  addition  of  a  higher  frequency,  even  thinner 
walled  pair  of  elements.  Lowering  the  driving  voltage  limit  has  the  effect  of  lengthening  the 
lower  frequency  elements  and  therefore  the  total  array.  Obviously  then,  shortening  the  full 
array  to  a  length  of  approximately  X/2,  as  recommended  for  the  prototype,  will  probably 
not  be  feasible. 

Another  approach,  however,  does  seem  feasible  for  controlling  the  directivity  charac¬ 
teristics  of  the  full  array.  The  requirements  for  the  flat  TVR  and  the  broad  main  beam  and 
low  sidelobes  are  contradictory;  that  is,  the  flat  response  requires  that  the  element  ampli¬ 
tudes  be  greatest  at  the  ends  of  the  array  and  decrease  toward  the  center,  the  beamwidth 
and  sidelobe  requirements  imply  just  the  opposite.  If  the  response  requirements  are  relaxed 
from  flat  to  a  positive  slope  as  a  function  of  increasing  frequency  (for  example,  +6  dB  per 
octave),  the  directivity  requirements  can  probably  be  met.  This  is  true  because  adding  a 
positive  slope  to  the  response  is  the  same  as  applying  a  linear  taper  (linear  shading)  to  the 
array.  The  output  sound  pressure  from  the  array  could  still  be  made  constant  as  a  function 
of  frequency  by  applying  the  reverse  slope  (-6  dB  per  octave)  to  the  amplitude  of  the  driv¬ 
ing  voltage.  The  computer  model  could  still  be  used  to  determine  the  element  sizes  required 
for  the  desired  slope  and  to  optimize  the  array  length  in  terms  of  the  directivity  and  effi¬ 
ciency  requirements.  The  positive  slope  approach  will  be  addressed  in  future  work. 
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